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<^ ; Abstract 

. The effect of the shock propagation on neutrino oscillation in supernova is studied paying 

,__( i attention to evolution of average energy of fe and Ue- We show that the effect appears as a 

CN ' decrease in average fe (in case of inverted mass hierarchy, Ue) energy at stellar surface as the 

^ ' shock propagates. It is found that the effect is significant 2 seconds after bounce if 3 X 10"^ < 

sin2%<10-2. 

^ : 1 Introduction 

^: 

^ Recently effects of shock propagation on neutrino oscillation in supernova was studied [H E] and it 
was shown that some characteristic signatures emerge as the shock propagates through the regions 
where matter-enhanced neutrino flavor conversion occurs. 



X 

^ ■ There have been many studies on neutrino oscillation in supernova: extracting information of 
neutrino parameters from the observation of SN1987A neutrinos |Sl IH EI or a future supernova 
neutrinos [3 IHl El HHl HI] , and probing supernova physics from observed neutrinos of a future super- 
nova But all of them are done without the effect of the shock propagation. 

In this paper the effect of the shock propagation is studied paying attention to evolution of average 
energies of z/g and Ue- We show when and with which parameter (sin^2^i3) the effect is significant or 
can be neglected safely. 

2 Neutrino oscillation and shock propagation 

If mixing angle is small, dynamics of flavor conversions is well described by resonant oscillation, 
which occurs at density, 

. , /Am2\ /10MeV\ /0.5\ 
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Figure 1: Evolution of density profile after bounce. Densities which correspond to H- and L-resonance 
are also shown. 

where Am^ is the mass squared difference, 9 is the mixing angle, Ei, is the neutrino energy, and 
Yf, is the mean number of electrons per baryon. Flavor conversion probabilities are determined by 
adiabaticity parameter 7: 

_ Aim? sin^ 29 
^ " 2E^ cos 2^ \dne/dr\' 

Here Am^ and 9 are 

^13 and Am^3 at H — resonance, 
9 12 and Am^2 at L — resonance, 

where mixing matrix is taken as: 

(C12C13 S12C13 Si3 \ 

-S12C23 - C12S23S13 C12C23 - S12S23S13 S23C13 , (3) 

■S12S23 ~ C12C23S13 — C12S23 — S12C23S13 C23C13 y 

where Sij = sin 9ij,Cij = cos 9ij for i,j = l,2,3(z < j). When 7^1, the resonance is nonadiabatic 
and the fluxes of the two involved mass eigenstates are completely exchanged. On the contrary, when 
7^1, the resonance is adiabatic and the conversion between mass eigenstates does not occur. 

Evolution of density profile is shown in Fig. ^ This is calculated by our numerical supernova 
model with progenitor mass, ISM©. (For detail, see JH]-) As the shock passes through the matter 
outside the neutrinosphere, densities will typically increase by approximately a factor of two. As 
neutrinos subsequently heat the matter behind the shock, a high entropy region will develop which 
accelerates matter behind the shock, causing a density buildup between this high entropy bubble 
and the outgoing shock. The density gradient caused by this mechanism is typically several orders 
of magnitude greater than that caused by the shock passage. As seen in Fig. ^ the shock is barely 
visible on the outer part of the density gradient at approximately 0.1 solar radius. This "snowplowed" 
density gradient is typical of neutrino driven core collapse supernovae, and is expected to be generic 
of a range of progenitor masses. 
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Densities which correspond to H- and L-resonance are also shown (Note that the resonance density 
depends on neutrino energy). As can be seen, the shock reaches the resonant region about 2 seconds 
after bounce and the density behind the shock is lower than that forward the shock. In the late phase 
[t > several seconds), a neutrino experience three times of H-resonances (and/or L-resonances) and 
the radii of resonance points (rres) are in general much smaller than the early phase. Since 7 oc r^es 
if the density profile is approximated to be power-law, in general the three resonances have different 
adiabaticity paremeters which are much smaller than that of the early phase. Therefore the average 
energy of the observed neutrinos is expected to depend on time not only due to the evolution of the 
neutrinosphere and the protoneutron star but also due to the evolution of the shock. 

We can estimate the above effect by calculating adiabaticity parameter. We consider the case 
where H-resonance occurs only one time and assume that the neutrino mass hierarchy is normal. In 
this case the adiabaticity parameter is, assuming p oc r"". 

Note that the index n is almost independent of time as can be seen in Fig. ^ From this H-resonace 
is expected to be adiabatic if sin^ 26^13 > 10~^. Adiabaticity parameter, however, becomes smaller as 
the shock propagates and two order smaller at 15 sec than in the early phase. Thus the H-resonance 
becomes less adiabatic unless sin^ 2^13 > 10^^ as is in the model LMA-L (sin^ 2^13 = 0.043) of 
[H]. This will cause decreasing of z/g average energy at the stellar surface as the shock propagates. 
On the other hand, if H-resonance is non-adiabatic even in early phase as is in the model LMA-S 
(sin^ 2^13 = 10~^) of FH|, average energy will not change. Here it should be noted that almost half 
the neutrinos are emitted after 2 seconds after bounce p!lj. 

With inverted mass hierarchy, H-resonance occurs at anti-neutrino sector. In this case evolution 
of average energy at stellar surface is seen in Ue- Its qualitative feature is expected to be the same 
as that of z/g with normal mass hierarchy because it is determined by the behavior of the shock. 

As to L-resonance, the adiabaticity parameter is, substituting the LMA (Large Mixing Angle) 
solution of the solar neutrino problem [121 CHI 113 UH] , 

Thus L-resonance will remain to be adiabatic even if the radius of the resonance point becomes two 
order smaller than the early phase. 

In the next section we study the above effect quantitatively by numerical calculation and obtain 
time evolution of the average energies of observed neutrinos (z/g and Ue)- 

3 Numerical calculation 

We solve numerically evolution equations of neutrino wave functions along the density profiles shown 
in Fig. nj From the wave functions, we obtain flavor conversion probabilities, from which neutrino 
spectra can be obtained by multiplying by the original neutrino fluxes. To make the shock prop- 
agation effect distinctive, the original energy spectra at each time are set to be the same as the 
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time-integrated spectra. Neutrino parameters are taken as: 

sin^ 2^12 = 0.87, sin^ 2^23 = 1, 

Am?2 = 7.0 X 10"^eV^ Am^g = 3.2 x lO^^eV^ (6) 

As for sin^ 2^13, we take various values including values corresponding to model LMA-L and LMA-S 
in IHlini- For detail of the calculational method and the original neutrino fluxes, see |Hlini- 

Fig. 121 show average energy evolutions of Vf. emitted at various times after bounce with normal 
mass hierarchy. Each figure differs in the value of sin^26'i3. The interesting behavior of the average 
energy in supernova (for example, 5 sec of the upper-left of Fig. El) indicates that H-resonance 
occurs three times. As is discussed in the previous section, the average energy at the stellar surface 
does not change in time when sin^ 2^13 is enough large (LMA-L, the upper-left of Fig. |21) or small 
(LMA-S, the lower-right of Fig. |21). On the other hand, average energy decrease by several MeV 
in the intermediate cases. By calculating with various values of sin^2^i3, we find that the shock 
propagation effects can be seen when 3 x 10~^ < sin^ 26*13 < 10~^ but are absent till ~ 1 second after 
bounce irrespective of sin^ 2^13. 

Thus average energy of observed Vf, will change in time due to shock propagating effect. In fact 
neutrino average energies changes also due to evolution of protoneutron star and neutrinosphere. 
Fig. 01 shows evolutions of average energy of observed neutrinos taking intrinsic changes of neutrino 
average energies into account, sin 2^13 is set to 10"^. As can also be seen in the lower-left of Fig. 
121 Vf. energy with shock effect is lower by several MeV after about 2 sec after bounce than without 
shock effect. 

As stated in the previous section, in case of inverted mass hierarchy it is z/g that is affected by 
shock propagation. Features of the evolution of the average energy are almost the same quantitatively: 
values of sin^ 26*13 and time after bounce, with which shock propagation effect is significant, difference 
between average energies of the early phase and late phase. 

4 Discussion and summary 

As we saw in the previous section, neutrino average energy decrease in general as the shock prop- 
agates. This is because the shock propagation cause decrease in the adiabaticity parameter of H- 
resonance and suppress the conversion between flavors. But it should be noted that the original 
spectra will change due to the evolution of the protoneutron star and the neutrinosphere. Thus we 
can not say about the value of sin^ 2^13 and the mass hierarchy only from the evolution of the average 
energy of the observed neutrino. To do so, we need a model of supernova and spectrum evolution. 

In this paper we studied the effect of the shock propagation paying attention to evolution of 
average energies of v^. and z/g. It is shown that the effect appears as a decrease in average Vf, (in case 
of inverted mass hierarchy, i/g) energy at stellar surface as the shock propagates. Further it is found 
that the effect is significant 2 seconds after bounce if 3 x 10~^ < sin^^is < 10~^. 



4 



Acknowledgements 

K.T.'s work is supported by Grant-in-Aid for JSPS Fellows. K.S.'s work is supported by Grant-in- 
Aid for Scientific Research (S) No. 14102004 and Grant-in-Aid for Scientific Research on Priority 
Areas No. 14079202. H.E.D. and J.R.W. are supported by LLNL contract W-7405-ENG-48. 

References 



[1] 

[2] 
[3] 
[4] 
[5] 
[6] 
[7 
[8 
[9 
[10 

[11 
[12 

[13 

[14 

[15 

[16 

[17 

[18 



R. C. Schirato and G. M. Fuller, |astro-p h/0205390| 



C. Lunardini and A. Yu. Smirnov, hep-ph/0302033| 



A. Yu. Smirnov, D. N. Spergel and J. N. Bahcall, Phys. Rev. D 49 (1994) 1389. 

B. Jegerlehner, F. Neubig, G. Raffelt, Phys. Rev. D 54 (1996) 1194. 
H. Minakata and H. Nunokawa, Phys.Lett. B 504 (2001) 301. 

C. Lunardini and A. Yu. Smirnov, Phys. Rev. D 63 (2001) 073009. 
A. S. Dighe and A. Yu. Smirnov, Phys. Rev. D 62 (2000), 033007. 

K. Takahashi, M. Watanabe, K. Sato and T. Totani, Phys. Rev. D 64 (2001) 093004. 
K. Takahashi and K. Sato, hep-ph/0205070 



K. Takahashi, M. Watanabe and K. Sato, Phys. Lett. B 510 (2001) 189. 

K. Takahashi and K. Sato, Phys. Rev. D, Phys. Rev. D, 66 (2002) 033006. 

H. Minakata, H. Nunokawa, R. Tomas and J. W. F. Valle, |hep-ph/01 12160) 

J. R. Wilson, R. Mayle, S. Woosley, T. Weaver, Ann. NY Acad. Sci. 470 (1986) 267. 

T. Totani, K. Sato, H. E. Dalhed and J. R. Wilson, Astrophys. J. 496 (1998) 216. 

S. Fukuda et al., Phys. Rev. Lett. 86 (2001) 5656. 

SNO collaboration, Phys. Rev. Lett. 89 (2002) 011301. 

SNO collaboration, Phys. Rev. Lett. 89 (2002) 011302. 

KamLAND collaboration, Phys. Rev. Lett. 90 (2003) 021802. 



5 




1-05 0.0001 0.001 0.01 0.1 1 10 

radius(solar radius) 



23 



CL> 



21 



2^19 
<D 

CD 
CD 

> 

CD 

13 
1e 



O.Ssec 


' ' 2sec 
5sec 

1; 15sec 

^1 1/ 


■ sin22ei3 = 10'^ 


/ )] 
/ /' f 

.--■/ 1 ^ L 







05 0.0001 0.001 0.01 0.1 1 10 

radius(solar radius) 




!-05 0.0001 0.001 0.01 0.1 1 

radius(solar radius) 



10 



> 



D3 
CD 

CD 

CD 
O) 
CD 

(D 
> 
CD 



19 

17 

15 

13 
1e 



sin220i3 = 10"^ 


...15sec 


/1 Osec '. 
.' 

; ■' \ 
,' / 

; / I 




• i 
// 

,v' .. — 


\ 
'\ 

5sec \ 




O.Ssec 



05 0.0001 0.001 0.01 0.1 1 

radius(solar radius) 



10 



Figure 2: Average energy evolutions of z/e emitted at various times after bounce with normal mass 
hierarchy. The values of sin^2^i3 are that of the model LMA-L of jHUH], 10~^, 10~^ and that of the 
model LMA-L of |Hl E] , respectively . 
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Figure 3: Evolutions of average energy of observed neutrinos. Average energies of Vf. with and 
without shock effect are shown. Those of z/g and — f^, i/^ without neutrino oscillation are also 
shown, sin^ 26*13 is set to 10~^. 
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